Aggregates of rod-coil diblock copolymers adsorbed at a surface 
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The behaviour of rod-coil diblock copolymers close to a surface is discussed by using extended 
scaling methods. The copolymers are immersed in selective solvent such that the rods are likely to 
aggregate to gain energy. The rods are assumed to align only parallel to each other, such that they 
gain a maximum energy by forming liquid crystalline structures. If an aggregate of these copolymers 
adsorbs with the rods parallel to the surface the rods shift with respect to each other to allow for the 
chains to gain entropy. It is shown that this shift decays with increasing distance from the surface. 
The profile of this decay away from the surface is calculated by minimisation of the total free energy 
of the system. The stability of such an adsorbed aggregate and other possible configurations are 
discussed as well. 

PACS numbers: 82.35.Gh Polymers on surfaces; 82.35.Jk Copolymers, phase transitions, structure; 36.20.Ey 
Conformation (statistics and dynamics) 
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I. INTRODUCTION 

Rod-coil copolymers in selective solvents show a rich 
phase behaviour. Depending on the chain length and 
the solvent quality they may form cylindrical micelles or 
lamellar sheets. Similar phases can be found in melts of 
rod-coil copolymers. These systems are therefore widely 
studied in the literature 0, H S S H B S S, H US . Less 
extensive are the studies on a single multiblock polymer 
composed of stiff rods which are connected by flexible 
chain spacers, see jO, 0, ^3 It was shown, that 
these polymers can also form micellar and multi-micellar 
structures. The structural behaviour of dissolved rod-coil 
copolymers in the presence of a surface is far less under- 
stood. Rod-coil polymers grafted to a repulsive surface 
are shown to form 'turnip'- or 'jellyfish'-like micelles on 
top of the surface However, to our knowledge there 
exists no study of rod-coil polymers in the presence of an 
attractive surface. 

In this paper we consider rod-coil diblock copolymers 
in selective solvent close to a surface which is highly at- 
tractive for the rods and neutral to the flexible parts of 
the copolymer. Further, the solvent is assumed to be 
poor for the rods, such that they align and tend to form 
aggregates, and good for the chains. In addition it is 
assumed that the rods have a certain chemical modifi- 
cation, such that they prefer to be parallel oriented with 
respect to each rather than antiparallel. The aggregation 
behaviour of such rod-coil copolymers, showing parallel 
alignment of the rods only, has been investigated exper- 
imentally and computationally, see 0, 0, |3 ■ We as- 
sume the energy penalty for antiparallel alignment of two 
rods to be much higher than the energy penalty for these 
rods being fully exposed to the solvent. In aggregates 
of these copolymers the flexible parts therefore stick out 
in one direction only, see Fig.Q. If such an aggregate 
adsorbs with the rods parallel to the surface, the rods 
shift with respect to each other to allow for the chains to 
gain entropy, see Fig.lQ). The nature of this shift will be 
examined in the following. 

For simplicity we consider a quasi two-dimensional sys- 



tem. That means the width of the system in y-direction 
is equal to the rod diameter d. This system can be viewed 
as a narrow slice of a system with infinite extension in 
y-direction. Each of the diblock copolymers under con- 
sideration is composed of a stiff rod of length L and diam- 
eter d to which a fully flexible chain of N monomers with 
monomer size h is grafted. The solvent is characterised 
by an energy penalty 7 per unit area of a rod exposed to 
the solvent. The energy gain —k per unit area of a rod 
for being in contact with the surface has to be chosen 
such that an aggregate actually adsorbs to the surface 
without dissociating into single rod-coil copolymers. 

The paper is organised as follows. In a first naive ap- 
proach one would assume a constant shift of the rods 
with respect to each other. This assumption leads to an 
artefact in the shifting behaviour, as we show in the ap- 
pendix. The calculations in the appendix are presented 
in some detail because some of the intermediate results 
are used in section II and III. In a more sophisticated ap- 
proach the shift is allowed to vary with distance from the 
surface such that it can develop a profile. This approach 
to the problem is presented in section II. Some remarks 
on the stability of the adsorbed structure and the corre- 
sponding range of values for k are made in section III. In 
section IV we finish with a brief discussion of the results 
of the foregoing sections. 



II. PROFILE OF THE SHIFT 

When an aggregate of rod-coil copolymers becomes ad- 
sorbed by a strongly attractive surface, an additional con- 
finement for the corona of the free chains may introduce 
new effects. One possibility is that the entropy penalty 
of the confinement of the corona chains due to the wedge 
defined by the surface and the rod aggregate prevents 
adsorption. Another extreme effect can be the destruc- 
tion of the aggregate. In between, the balance between 
entropy and energy can be such, that the rods are shifted 
with respect to each other to allow for the chains to gain 
entropy. If this shift does not get to large, the rods stay 
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FIG. 1: Aggregate of rod-coil copolymers adsorbed at a sur- 
face. The shift of the rods with respect to each other decays 
with increasing distance from the surface. 



together. Such conformations appear possible whenever 
the defect energy is balanced by the gain of entropy due 
to reduction of the wedge confinement. 

The most naive assumption that the shift may be con- 
stant, i.e., it forms a line with a certain slope, turns out to 
be unphysical in many respects (see Appendix). There- 
fore it is necessary to introduce a curved deviation of the 
shift for each successive rod, which results from a total 
balance of all entropic and enthalpic contributions. 

We present in this section the basic model which al- 
lows us to calculate the shift of the rods as a function 
of distance from the surface. From intuition we expect 
the rods close to the surface to shift more than the rods 
further away from the surface, since for the correspond- 
ing chains close to the surface there is more entropy to 
gain than for the ones further away. Therefore we expect 
an equilibrium conformation similar to the one shown in 
Fig.©. 

The dashed line in Fig.Q can be interpreted as the 
profile of the shift I as a function of the distance x from 
the surface. The chains are described by a local Flory- 
type model similar to the one used in [13 to describe a 
finite polymer brush. In this model the free energy of 
the system is given by the sum of an elastic term and an 
excluded volume term for the chains plus a term which 
quantifies the energy penalty for the additional rod sur- 
face exposed to the solvent due to the shift. We construct 
the free energy such that it is a function of the splay of 
the chains u{x) - see Fig.Q - and the shift of the rods 
l{x). At all positions x we can safely assume l(x) to be 
small compared to the height h of the brush-like struc- 
ture formed by the chains. If the shift I would be of the 
same order as h, the chains would hardly see each other, 
so there would be no driving force for a shift. Hence 
we can assume h to be constant for all x. The height 
is given by the equilibrium height of a polymer brush in 
good solvent, i.e. h = A^''-/^b N{b/d)'^/^ . The excluded 
volume parameter v is set to v — for simplicity. 

Fig.||2Jl helps to define and explain how the excluded 
volume term of the free energy for the shift geometry is 
constructed. Note, that it is only a sketch. The ratio of 



shift I to brush height is much smaller than depicted in 
Fig-©- A chain starts at the rod and ends at the line 



X(x) 




FIG. 2: This sketch shows the spatial segments or boxes filled 
by each chain. They get larger with increasing distance from 
the surface. The shift at each position x is denoted by l{x). 
X ranges from —D at the surface to at the last rod. The 
splay of the chains is given by u{x) = X[x) — x. 

X{x) shown in Fig.®. It is assumed to fill the volume 
of the box given by the dashed lines around the chains. 
We are aware that this assumption is not valid for the 
chains far away from the surface but for these chains 
the contribution of the excluded volume term is certainly 
very small. Hence the assumption does not affect the 
total free energy in a significant way. The splay u{x) 
is given by X(x) — x. To explain how to calculate the 
volume available to each chain, Fig.(|3Jl shows a larger 
sketch of one of the dashed boxes surrounding each chain 
in Fig.©. 

Inasmuch as we are going to consider a quasi two- 
dimensional system only, the volume is given by the grey 
area in Fig.lO times the diameter of the rods d. The 
area A-°^ shaded in light grey is given by h{d+ Aui/2), 

where Au^ — Ui — Ui-i. It is the area A^"^ shaded in dark 
grey where the shift of the rods l{x) comes into play. It 
is given by 
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As can be seen from Fig.© the length qi is given by 
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FIG. 3; This figure shows how the volume occupied by each 
chain is calculated as a function of splay u and shift I. 



Since we assumed U to be small compared to h, the term 
in the last brackets in Eq.lPJ can be approximated by 1. 
In Fig.|(2Jl this corresponds to double counting the area 
of the small triangle with the catheti qi and U. 

In order to construct a free energy functional that 
could be minimised with respect to the shift and splay 
shapes l{x) and u{x) we take the continuum limit: li — > 
l(x), Ui u{x), Aui du'{x). Taking the continuum 
limit and adding the two contributions, the total volume 
V{x) available to a chain at position x is given by 

V{x) ^(fh{l + u'{x)/2)+dl{x)(u(x)~du'{x)). (4) 

The elastic term of the free energy is given by a contri- 
bution proportional to representing the stretching of 
the chains away from the rods and a contribution propor- 
tional to u{xY representing the stretching of the chains 
parallel to the rods. The energy penalty for the addi- 
tional area of a rod at position i exposed to the solvent 
as a function of the shift U is simply given by 2^lid. The 
complete free energy functional can now be constructed. 



-D 



2l3 



2d 



PF = 2 i dxPjl{x) + 



dx 



Nb^d 



dx [u{3 



\2 , l21 



-D 



hd^{l + u'{x)/2) -f dl{x){u{x) - du'{x)) 



The excluded volume term in Eq.lO was constructed 
assuming constant density of monomers for each chain 
within the box of volume V (Eq.Q). The monomer 
density close to the rods is larger than the one further 
away from the rods and therefore this assumption tends 
to underestimate the excluded volume energy. However, 
it is the standard approximation used in Flory-type mod- 
els and has been proven to be sufficient to describe the 
behaviour of a finite polymer brush, see . 

The total length D of the rod aggregate perpendicular 
to the surface is assumed to be D > h. This is explicitly 
needed in the appendix. 

Fig. 10) shows that the rod at the surface has zero shift 
since there is no other rod underneath with respect to 
which it could shift. So for the rod-coil copolymer at 
the surface the integrand in Eq.jSJ reduces to the one in 
Eg. ljAip in the appendix (with a = l/rf^). This of course 
also means that the equation for the shift (Eq. (|SJ| below) 
is only valid from the second rod on (as counted from the 
surface) . 

To calculate the equilibrium shift l{x) it is necessary to 
compute the Euler-Lagrange equations from a functional 
minimisation of Eq.|SJl with respect to u{x) and l{x). The 
Euler-Lagrange equation for the splay u{x) has the first 
integral 
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u{xY + h"^ 
Nb'^d 



hd^{2 + u'{x)) + 2dl{x){u{x) - du'{x)) 
{h-2l{x))u'{x) 



{hd{2 + u'{x)) + 2l{x)(u{x) - du'{x))Y 



= C,. (6) 



Variation of the free energy functional in Eq.(|S)) with 
respect to l{x) yields 



u{x) — du'[x) 



{hd(2 + u'{x)) + 2l{x){u{x) - du'ix)))' 
The quadratic Eq.© can be solved for the shift l{x). 



(7) 
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(8) 



For simplicity we dropped all terms which do not depend 
on either I or u. As already mentioned the excluded 
volume parameter v is set to v = b^. 



Inserting the expression for the shift - i.e. Eq.® - into 
Eq.® results in a complicated, highly nonlinear differen- 
tial equation which cannot be solved exactly. However, 
the overall effect of the shift on the free energy is cer- 
tainly smaller than the overall effect of the splay. Thus 
it is a reasonable approximation to calculate the solution 
for the splay at zero shift and to use this as an approxi- 
mation for the splay u in Eq. (jHJ . We are aware that this 
(5) approximation breaks down when the splay becomes very 
small close to the surface. Nevertheless, for zero splay the 
shift should be constant. The shape of the profile of the 
shift away from the surface at finite splay can therefore 
be calculated within this approximation. 
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Now we calculate an approximate solution of the splay 
u for zero I. If the shift l{x) in Eq.lQ is set identical zero, 
the differential equation equation reduces to 



2N^ 



Nb^d 



l + u'{x) 



(2 + 
8/1^ 



u'{x)y 
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(9) 
(10) 



which is Ea. l|A2(l from the appendix as should be. It 
is convenient to introduce dimensionless variables u = 
m/(V8/i), X = x/{V%h) and C2 = C2/(8/i2). Eq.lTOI) can 
be integrated for arbitrary C2 which leads to an implicit 
equation for the splay u{x) similar to Ea. HA3|l 



2uq - 2w + In 2m + (1 - 4(72 + ^u'^f'^ 
- In [2uo + (1 - 4(72 + 4ug)^/2j = 4i 



(11) 



Note that —D < x <{). The integration constant C2 can 
be determined by using the boundary condition u[—D) — 
- i.e. zero splay at the surface, see Fig.(|21l. 
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(12) 

The chain furthest away from the surface at a; = can 
topple over completely and is therefore allowed a splay 
Uq = h OT Uq = l/\/8. Therewith Eas. (|llll2|l reduce to 



l/\/2- 2{t + ln 



2u + (1 - 4(72 + 4{t^) 



^,2^l/2 



-In 
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(72 = - I 2 - sinh 



AD + l 



= 4i, 

-2^ 



(13) 
(14) 



It is not possible at this stage to resolve Eq. (|13|l with re- 
spect to the splay u. However, we know function 
of u, u' as a function of u (see Ea. (|10|) ') and I as function 
of u and u' . Hence for a certain set of parameters, we 
can plot the shift I as a function of x by either numer- 
ically resolving Ea. (|13|l or by showing a parameter plot 
of I versus x, using m as a parameter. This is done for 
a characteristic set of parameters in Fig.l^J. The plot 
demonstrates that our assumption for the profile of the 
splay - Fig.J^I - was indeed reasonable. The shift in- 
creases from the furthest rod towards the surface until it 
reaches its maximum value in a reasonable form. 

We are now going to estimate the threshold value of 7 
above which the energy penalty for additional rod-solvent 
exposure becomes to large for a shift to occur. 

The shift is identical zero if the right hand side of 
Eq.® is less or equal to zero for all x. 
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FIG. 4: The shift I is plotted as a function of negative dis- 
tance. The rod furthest from the surface is located at a; = 0, 
compare Fig.lUJ. Parameters: iV = 700, b/d = 0.1, /S-fd"^ = 1 



To find an upper limit for the threshold value 7c we note 
that the maximum value of the shift is uq = h. An upper 
estimate for {u{x) — du'{x))/{2 + u'{x))'^ is therefore given 
by Ah. The upper limit for 7c is thus given by 



/37c 



b^N^ 
Ahd^ 



0.4 



.N 

d2 



(16) 



For the set of parameters chosen in Fig. this yields the 
rough estimate of [3"fcd^ « 13. 

In this section an attempt was made to describe a 
variable shift l{x). For values of the splay u{x) which 
are large enough to dominate the effect of the shift, we 
found a set of equations JH] El I14|l which determine I as a 
function of u and u as a function of x. Although it is not 
possible to resolve Ea. H13|) with respect to w, the profile 
of the shift can be plotted for a certain set of parame- 
ters, see Fig.(0J. In this section we always assumed that 
the adsorbed aggregate is internally stable and does not 
disintegrate in the sense that single copolymers leave the 
aggregate. Therefore we are going to discuss in the next 
section under which conditions the adsorbed aggregate is 
stable and which other configurations are possible. 



III. STABILITY 

In section II we assumed that an attached configura- 
tion of / rod-coil copolymers at a surface as it is pictured 
in Fig.(^) is stable. There are three other possible con- 
figurations. 

In this quasi two-dimensional system only one rod is 
in contact with the surface. Therefore a possible con- 
figuration is the one shown in Fig.©, where one single 
copolymer is adsorbed at the surface and the others form 
a detached sheet. We call such a situation detached con- 
figuration in the following. The rods in the detached 
sheet might also prefer a shifted geometry. We refrain 



from a discussion of this shift, since this paper is mainly 
concerned with the behaviour of a complete aggregate in 
contact with a surface. 



FIG. 7: The aggregate dissociates and the individual rods 
adsorb at the surface. This configuration is preferred for large 




FIG. 5: For long chains a detachment of the aggregate from 
the rod adsorbed to the surface might be preferable. 

The energy of the rod-surface contact in the detached 
configuration is the same as for the attached one. A con- 
figuration with different contact energy is the mushroom- 
like one as depicted in Fig.JHll. This configuration is al- 




small the system might prefer the mushroom configura- 
tion even for long rods, since it allows for the chains to 
gain entropy without much increase in exposure of the 
rods to the solvent. 

By estimating the free energies of these configurations 
and comparing them with the one of the attached con- 
figuration, it is possible to find the range of k in which 
the attached configuration is stable. To achieve this at 
least approximatively we calculate the free energy of the 
attached configuration with zero shift, ft gives a slight 
overestimation of the free energy of the configuration 
considered in section 11. However, we still get a rough 
estimate for the parameter range in which the attached 
configuration is stable. At zero shift the chains form a 
finite brush. Its free energy is calculated in the appendix 
and given by Eq. (jAll|l . One side of the brush is free and 
therefore allowed a splay of uq = the other side is con- 
fined by the surface, i.e. ub — 0. The length D is given 
by fd. The free energy of the attached configuration then 
reads 
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attach 



Nb^d 



3fdh^ ^ {- - ^] 



(17) 



We use this as a reference energy and add energy gains 
and penalties due to rod-surface contact or rod-solvent 
exposure to the free energies of the other configurations. 

Compared to the attached configuration the contact 
energy of the mushroom differs by K{Ld — fd^). The 
chains can also be described as a finite brush. Here both 
ends are free and are therefore allowed a splay of uq = 
ud = h. The free energy of the mushroom configuration 
is hence given by 



FIG. 6: Another possible configuration: The aggregated rods 
adsorb perpendicular to the surface. This configuration is 
preferable for large aggregates. 

ways preferable to a complete detachment of the aggre- 
gate since in the latter case the system would gain no 
contact energy. 

The last possible configuration is a complete dissocia- 
tion of the aggregate into single copolymers due to the 
presence of the attractive surface. These single copoly- 
mers then adsorb individually at the surface, see Fig.lQ. 
This configuration yields the highest gain of contact en- 
ergy. However, it is also the configuration with the high- 
est energy penalty for exposure of rod surface to the sol- 
vent. For very high contact energy, i.e. k ^ 7, the 
system always dissociates. On the other hand, if n is to 



musnroom 



K (dL - fd^ 



1 



Nb'^d 



■ifdh^ 



.(18) 



The attached configuration is preferred to the mushroom 
configuration if Fattach < Fmushroom- This yields the 
following condition for k 



/3k > 



0.127V2 

Ld - fd^ 



(19) 



This is the lower bound for k. To get the upper bound 
the free energy of the dissociated copolymers - see Fig.© 
- has to be estimated. 

Compared to the attached configuration the dissoci- 
ated one yields a contact energy gain of — (/ — 1)kLo?. 
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But on the other hand it also gives rise to an additional 
energy penalty of 2(/ — l)^Ld. Within this Flory-type 
theory the flexible chains of the individual copolymers at 
the surface can be treated as free ones and their free 
energy can be neglected. Comparison of Fattach and 
2(/ — l)jLd — (/ — l)KLd yields the upper limit of k 
above which the system dissociates: 



detached configuration has a contribution from two addi- 
tional rod surfaces exposed to the solvent, which is given 
by 2Ldj. The free energy of the chains of the detached 
sheet is similar to the one of the mushroom configura- 
tion, with D = / — 1. The free energy of the chain of 
the single copolymer can be neglected as in the case of 
dissociation. Fdetach is hence given by 



(5k < 2/37- 



(20) 

Within this range of values for k the attached configura- 
tion can actually be stable. For the parameters chosen 
in Fig.Q, N = 700, b/d = 0.1, (jjd"^ = 1 and L/d = 80, 
/ = 30, this range is given by 1.18 < /3kc?^ < 1.52. 

It is also of interest to keep k and 7 fixed and to in- 
vestigate at which combination of molecular properties 
of the copolymers which configuration is preferred. We 
focus here on the length of the rods L and the number 
of chain monomers N. In the following the critical rod 
lengths which separate each two of the possible config- 
urations from each other are calculated as functions of 
N and of the other parameters. Since it depends on the 
specific combination of parameters, especially on 7 and 
K, which configurations are neighbouring in L-N space, 
all possible critical rod lengths are calculated. Note, that 
not all possible phase boundaries can exist for one given 
set of parameters. For each one to exist, the set of pa- 
rameters have to be chosen appropriately. 

Ea. H19|l can be rearranged such that it gives the criti- 
cal rod length below which the system changes form the 
attached to the mushroom configuration. 



T attach 
mushroo7n 



fd + 0.12 



(3Kd^ 



(21) 



Comparison of Fj 



^mushroom with the energy of the dissoci- 
ated configuration yields the rod length L'^^^^J^""^ which 
separates the mushroom from the dissociated configura- 
tion. 



L 



musnroora 
dissociate 



1.2/diV(^)'/' 



0.2AN%{^f ~ fd[iKd^ 



2{f - l)pjd^ - ff3Kd^ 



(22) 

In case of rather large k (close to the upper limit) there 
might exist a rod length which directly separates the at- 
tached configuration from the dissociated configuration. 
It is found to be 



J- attach 
dissociate 



l.2fdN{^,f'- 



0.12N%{^y 



if - l)(2/37d2 - pud^) 



(23) 



To calculate the boundaries of the detached configura- 
tion (see Fig.©), its free energy has to be estimated. As 
already mentioned above, this configuration might also 
prefer a shifted geometry. But since we already esti- 
mated Fattach for zero shift, is is sufficient to calculate 
also Fdetach for a rectangular sheet without shift. Com- 
pared to the attached configuration the free energy of the 



(iFdetach = 2Ldp-i+ 



1 



Nb^d 



3(/ - l)dh' 



(24) 

The rod length L'j^uach with separates the detached and 
the attached configuration can now be estimated. 



T detach ^ n 
^attach ~ "-O 



N 



df3j \d 



d(3j \d 



However, in L-N space the detached configuration might 
sit in between the mushroom and the dissociated config- 
uration. This is indeed the case for a wide range of pa- 
rameters. Therefore equating Ea. H18l) and Ea. H24l) gives 
the rod length which separates mushroom and detached 
configuration. 



- musriroom 
^detach 



fd- 1.2 



Nd 

l3Kd^ 



1-1 



(26) 



The length L'^fg^o^^^te which separates the detached from 
the dissociated configuration is found to be 



1.2(/-l)cijV(^)'/'-0.24jV^b(^) 
2(/ - 2)/?7d2 - (/ - l)(3Kd') 



T detach 
dissociate 



(27) 

In the conclusions we use these critical rod lengths to 
calculate one phase diagram in L-N space for a typical 
set of parameters as an example. 

In this section it was shown that the configuration pic- 
tured in Fig.J^I can be stable within a certain range of 
values for k. There are three other possible configura- 
tions as shown in Fig.©, Fig.© and Fig.©. 



IV. CONCLUSION 

A discussion of rod-coil copolymer aggregates adsorbed 
at a surface in a two dimensional approximation was pre- 
sented. The aggregates form because of a selective sol- 
vent, poor for the rods and good for the chains. Due to 
their chemical structure the rods only align parallel to 
each other. The surface is assumed to be attractive for 
the rods and neutral with respect to the chains. If the 
aggregate adsorbs with the rods parallel to the surface, 
the rods shift with respect to each other to allow for the 
chains to gain entropy and to therefore lower their con- 
finement energy. If the shift is assumed to be constant for 
all rods, the system shows an artificial behaviour. This 
can be seen from the considerations in the appendix. 
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In section II we constructed a model which allows for 
the shift to vary. It was possible to partially solve this 
model and to show that away from the surface the shift 
decays. Close to the surface the splay of the chains is 
close to zero and therefore the shift is basically constant. 
The region further away from the surface is the interest- 
ing one showing the decay profile of the shift, see Fig.Q). 

In section III we showed that the configuration con- 
sidered in section II can actually be stable within a cer- 
tain range of contact energies between rods and surface. 
This range was calculated. Three other possible con- 
figurations were also discussed, see Figs. (|5l6l7|l . These 
considerations allow us to plot a phase diagram of the 
configurations in L-N space, which - for a typical set of 
parameters - is shown in Fig.JHJ. The contact energy 
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FIG. 8: Configurations of the rod-coil copolymers at the sur- 
face. Parameters: b/d = 0.1, / = 30,/37d^ = l,/3Kd^ = 1.3 

per unit area k is chosen such that there exists a region 
in L-N space in which the entropy loss of the confined 
chains is compensated by the energy gain due to rod- 
surface contact. But it is also chosen to be not much 
larger than the rod- rod contact energy 7, since other- 
wise the aggregate would dissociate. However, for very 
long chains the aggregate always dissociates into single 
copolymers which then individually adsorb. Neverthe- 
less, Fig.® shows that there is indeed a broad region in 
L-N space in which the attached configuration is stable 
and the rods shift with respect to each as discussed in 
section II. 

For further considerations in the future it would be of 
interest to study the behaviour of a finite three dimen- 
sional aggregate adsorbed with the rods parallel to the 
surface. We expect the profile of the shift to form a two 
dimensional surface with the innermost rods close to the 
surface showing the maximum shift. 



APPENDIX A: CONSTANT SHIFT 

Despite the fact that the assumption of a constant shift 
appears unphysical and leads to contradictory results we 
are discussing it in this appendix in some detail. It is a 



useful example to define the arising problems in a clear 
way. 

The rods are assumed to shift a constant distance with 
respect to each other as shown in Fig.©. The charac- 
teristic quantity related to this shift is the angle a. This 
angle can be calculated by calculating the free energy of 
the entire system and minimising it with respect to a. 

The additional free energy per rod due to the shift is 
given by the additional surface of the rod exposed to the 
solvent Frod = 27d^tana. To calculate the free energy 
of the chains we treat them as if they would form a fi- 
nite brush grafted to the surface shown as a thick line in 
Fig.®. For a finite brush the trajectories of the single 
polymer chains are not all perpendicular to the grafting 
surface as for an infinite one. The polymer chains show a 
splay u. Fig.® illustrates how this quantity is defined; 
that is, u{x) = X{x) — x. The first chain is allowed a 
splay of uq — htana due to the surface, where h is the 
brush height given by ft. = A~^^^b Na^^^b'^^^ . The last 
chain can topple over completely and is therefore allowed 
a splay ul = h. The grafting density tr is a function of a 
as well. It is given by ij{a) = cos{a)/cP. As in section II 
a Flory-type approach is used to describe the free energy 
of the finite brush following the lines of . Each chain 
fills a box of volume ha^^{l + u' /2), with u' = du/dx. 
The free energy for the finite brush is then given by 



brush 



Nb"^ 



dx [u{x)^ -\- h 



D 






1 




1 dx 






2 + u\ 





(Al) 

where D represents the total length of the brush. For 
our system of aggregates of rod-coil copolymers it is 
given by D{a) = /dcos(a). The first integral of the 
Euler-Lagrange equation for the splay u{x) obtained from 
Eq. (|Aip is given by 



8/1^ 



1 



C. 



(A2) 



(2 + u')2 

Eq. (|A1|I does not distinguish between positive and nega- 
tive splay. Hence we have to separate our system into two 
finite brushes which meet at the chain with zero splay, 
see Fig.®. This chain has to be determined by an equi- 
librium condition (total length of the brush D). 

It is convenient to introduce dimensionless variables 
u = u/{VSh),i; = x/iVSh) and C = C/(8ft2). Ea.(IX2ll 
can be integrated for arbitrary C which leads to 



I = 2-u - 2uD - In 2-it + (1 - 4C + 4u 



In 



2uD + {I - 4,0 + Auiy/^ =A{b-io-i). 

(A3) 

The integration constant can now be calculated by using 
the condition 2I(u 0) — l{u = uq) 

1 



(A4) 



C = -l^l-iuo+uo) sinh[2D + uo + UDY 
-h {uq - ud)"^ cosh[2D + Uq + udY 
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then be rewritten in the following form 




FIG. 9; Aggregate of rod-coil copolymers adsorbed at a sur- 
face. The rods are all shifted with respect to each other by the 
same distance. The shift is characterised by the angle a. This 
sketch shows how we define the a;-range and the splay. Note: 
u{x) — X(x) — X, i.e. uq = Xo — xq and ud = Xd — {D — xq). 



The free energy in Ea. ljAl|l cannot be integrated di- 
rectly using the impHcit solution for the splay u, Ea. (jA3p . 
Therefore we have to find an appropriate approximation. 
It can be shown that 1 — 4C w for all uq ii D > h. 
Hence Ea. l|X3|) can be very well approximated as 



2u — 2uo — In 
for 



UD 



4:{D-io - i), 



u>Uc^l{l~4Cy/^. (A5) 



For < u < Uc we can choose a linear approximation. 

To calculate the free energy we have to perform the 
following integration 



da; 



1 



2(2 + u') 

We first consider the regime < u < Wc 



u{x < Xc) < 2 



1 



2(2 -I- u') 



1 

4' 



(A6) 



The integral in Ea. l|A6p in the interval [0, ic] can there- 
fore safely be approximated by 



dx 



1 



2(2 + u') 



1 ^ 

-Xc 



(A8) 



In the interval [x^, D — xq] the approximation for the 
splay, Ea. (|A5p . is valid. The integral in Ea. ljA6|) can 



D-Xa 

I 



dx 



1 



2(2 -t- ii') 



dx 



uu' — u^u' 



u' 1 
4+2 



X u 
4 ~ 8 



~2 -■^T-D- 



6 



(A9) 



By construction u{xc) = Uc ^ 1, see Ea. (|A5|) . Therefore 
the integral in Ea. (jA6|) in the limits [0,13 — io] can be 
very well approximated by making use of Eas. (|A8IA9p : 



D-5 



dx 



I? 



2(2 + u') 



D 



Xq Ud 



4 



6 ' 



(AlO) 

So far we calculated only one part of the brush. The 
one from the "zero splay chain" to the open end. In 
Fig.® this is the right part from to D — xq- The left 
part from to xq or rather from the "zero splay chain" 
to the surface can be calculated completely analogous 
replacing D — xq with xq and ud with uq. Adding up the 
results for both parts in both intervals, accounting for the 
prefactors in Ea. l)Al(l and converting back to variables 
carrying dimensions we get as a final result for the free 
energy of the chains forming the finite brush 



brush 



a \iDh^ - h^{uo + UD) 



h ( 2 



(All) 



Plugging in the a-dependent expressions for cr, D, uq and 
Ud and adding Fj-od we get the a-dependent part of the 
total free energy of the system. 

/3F(a) = 2f(3jd^ tan(a) -f -^fN ( - \ cos(a)^/3 



cos^(a) X 



71"^'*'"^"^^ V2 



4 U. 

tan^(Q!) tan'^(Q;) 



6 



(A12) 



This equation is only valid for < a < 7r/4, since 
uo(7r/4) = /i = In the interval 7r/4 < a < tt/2 

the splay at the surface remains constant at its maxi- 
mum value Uq — h, and the term in square brackets in 
Ea. (|Il2)l reduces to [\/2 - 7/3]. 

There are two regimes in each of which the above 
free energy, Ea. (|A12|l . shows a different behaviour. The 
regime where the chains are rather short shows a first 
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order like transition. With decreasing 7 there is a jump 
from a stable phase with no shift {a = 0) to a stable 
phase with a large shift (a 3> 0). The other regime in 
which the chains are rather long shows a second order 
transition ft'om a stable phase with a = to a phase 
with finite a. 

As a criterion to distinguished these two regimes the 
chain length can be used. The chain length Nt which sep- 
arates these two regimes is given by Nf w 1.4/ (f )^^'^- 
For N < Nt the transition is first order like and con- 
tinuous for N > Nf. This condition can be interpreted 
such that for N < Nt the second term in the free energy 
fEa. (|A12|l ^ dominates the third term. The second term 
solely represents the effect of decreasing grafting density, 
whereas the third term also represents the effect of in- 
creasing splay of the chains close to surface. Therefore, 
if the grafting density effect is dominant, the transition 
is similar to the tilting transition observed in lamellar 
structures of rod-coil copolymers, see e.g. 0. Since then 
the minimum in the free energy is essentially given by 
the balance of tan a (first term) and cos a (second term) 
there cannot be a continuous transition with a minimum 
at small values of a. To illustrate this behaviour, the 
free energy as a function of a for different values of 7 is 
plotted in Fig. lfTUjl . 



F 




40 - ^ _ / 



' ' ' — — ^ — — — a 

0.2 0.4 0.6 0.8 1 1.2 1.4 

FIG. 10: This plot shows the dependence of the free energy in 
units of kT on a for different values of 7 in the regime N < Nt- 
Parameters: / = W,b/d = 0.2, AT = 50, /3^d^ = 1 (upper 
curve), 0.25, 0.15 (lower curve). In this regime the effect of 
decreasing grafting density with increasing a dominates the 
effect of the splay of the coils close to the surface. 



For N > Nt the third term in the free energy, 
Ea. (|A12|l . gets equal to or bigger than the second term. 
This means that the increase in splay of the chains close 
to the surface becomes important. Since uq scales with 
tana like the contribution of the rods (Frod) does, a con- 
tinuous transition is now possible. The dependence of 
the free energy on a for different values of 7 is shown in 

The assumption of the shift to be constant is an over- 
simplification. The crossover from one regime {N > Nt) 
in which a shift develops continuously to a regime {N < 
Nt) in which the system shows a first order transition like 
jump from zero shift to large finite shift is an unphysi- 
cal artefact of this approximation. Only in the limit of 
a very large number of copolymers / forming one lamel- 
lar like aggregate the constant shift assumption might be 
reasonable. However, in this limit the effect of the sur- 
face becomes negligible and the system always shows a 
tilting transition - compare |^ - even if not in contact 
with the surface. 



F 




a 



0.2 0.4 0.6 0.8 



FIG. 11: This plot shows the dependence of the free energy 
in units of kT on a for different values of 7 in the regime 
N > Nt. Parameters: / = 10, b/d = 0.2, = 200, f3-yd^ = 4 
(upper curve), 3, 2.5 (lower curve). In this regime the effect of 
increasing splay of the chains close to the surface determines 
the behaviour of the free energy. 
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